NIKHEF-94-P9 
NPL-1110 



|hep-ph/ 941 2287 



Evolution of the spin of the nucleon 

P.J. Mulders * 

National Institute for Nuclear Physics and High Energy Physics (NIKHEF). 
P.O. Box 41882, NL-1009 DB Amsterdam, The Netherlands 
Email address: pietm@nikhef.nl 

S.J. Pollock 

University of Colorado, Boulder CO, USA 80303 
Email address: pollock@lucky.colorado.edu 

Abstract 

We compare momentum sum rules from unpolarized electroproduction and 
the spin sum rule for gi in polarized electroproduction, and their evolution 
in the framework of the operator product expansion. Second order effects in 
as are included. We show that in comparing the evolution of the spin sum 
rule with the momentum sum rule one is not overly sensitive to using first or 
second order, even when going to the extreme low limit in which gluons 
carry no momentum. Our results show that in that limit there is no need to 
include any contribution of strange quarks. 



*Also at the Department of Physics and Astronomy, Free University, Amsterdam 

1 



I. INTRODUCTION 



Deep inelastic scattering (DIS) is an important tool for studying the structure of liadrons. 
Through sum rules the experiments provide values of specific quark and gluon operator 
matrix elements. The framework for this is the operator product expansion (OPE) 
Experimentally measurable sum rules are expressed as the product of matrix elements and 
coefficient functions. Examples are the momentum sum rules measured in unpolarized deep 
inelastic scattering and the Bjorken sum rule and Ellis- Jaffe sum rule |^ in polarized 
deep inelastic scattering [Q. Initial measurements of the latter, showing deviations from 
the Ellis- Jaffe prediction, have been interpreted as an indication of a surprisingly small 
contribution of the quark spin to the nucleon spin One of the points relevant for the 
interpretation is the scale dependence of the matrix elements and the coefficient functions, 
which can be calculated in perturbation theory. The QCD corrections to the Bjorken sum 
rule up to order have now been calculated in leading twist and the higher twist 

corrections have been estimated 0. Recently, the order corrections to the Ellis- Jaffe 
sum rule in leading twist and massless quark approximation have been completed as well 
0. These corrections provide powerful means to further study the evolution of the spin 
structure of the nucleon. 

For matrix elements that have no scale dependence deep inelastic measurements imme- 
diately provide us with interpretable results that occasionally can be compared with other 
experimental measurements in a completely different domain, e.g. the Bjorken sum rule. It 
is well known ^ that the singlet part of the first moment of the spin structure function gi 
is not of this type and has an anomalous evolution. It is also well known |JlO[ that the 
leading term in as in the axial anomalous dimension vanishes, and for this reason some au- 
thors dismiss this evolution as insignificant. Roughly speaking, corrections to the singlet 
first moment arising from the anomalous dimension can be argued to behave like as logQ^, 
and hence appear approximately independent. In an earlier paper, |]ll| we argued that 
comparing momentum sum rules from unpolarized electroproduction and the spin sum rule 
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for gi, including their evolution, showed that DIS spin measurements are consistent with 
a low energy valence quark picture, where the valence quarks carry a substantial part of the 
spin of the proton, namely of the order of Ga/{5/3) ~ 0.75. In this paper, we extend our 
earlier calculations to fully include the next higher order QCD corrections in leading twist. 
In this way we are able to get a feeling for the sensitivity to the use of first and second order 
in the evolution of the spin sum rule. We also consider the possible effects of a strange quark 
threshold at very low momentum scales. With the results of more recent experiments we 
can give an estimate of the "spin carried by quarks" . 



II. FORMALISM 



A. The momentum sum rules 

As a typical example of an (unpolarized) "quark momentum sum rule", the second 
moment of F2 is given by 

x[qi{x, Q"^) + qi{x, Q"^)] dx 



I dxF2ix,Q')=J2eU 
Jo J 

= T.ehm = Ee^r(Q^) + {e')MQ% (1) 



i=l i=l 



where qi{x, Q ) is the quark distribution function and ej = ei{Q ) is the momentum fraction 
carried by quarks and antiquarks of flavor i (n^ is the number of flavors), which is separated 
into nonsinglet (NS) and singlet contributions. The quantity S2 = J2i U is the total mo- 
mentum fraction carried by the (valence+sea) quarks, which can be expressed as the matrix 



element of the quark part of the energy momentum tensor. The quantity ef = e.j — ^2/^/ i 



IS 



the nonsinglet part of the second moment for flavor i, and (e^) is the average quark charge. 
There is no unique way to define parton distributions beyond leading order, but we follow 



Buras' |jT2[ conventions, including renormalization in the MS scheme. This results in the 



following formulae for the (unpolarized) momentum sum rules, including next to leading 



order corrections ]T3[. For the (nonsinglet) valence quark momentum sum rule 
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^2 = E /' x[q,{x, Q^) - q^{x, g2)] dx (2) 
the evolution is given by 



with the anomalous dimension given by 



7™(".)=7.-(^)+ -,f' (4) 

with 7,^'^ = 64/9 and 'y^^ = 96.6584 — 6.32 n/. The beta function governs the behavior of 
the strong coupling constant, 

,^0=,(„.)._,„g_ft_^_..., (5) 

with (3o = 11 - 2n//3, (3i = 102 - 38n//3. 

The leading order (LO) solution for the strong coupling constant is 

^^^-log(g^)= invariant, (6) 

while the next to leading order (NLO) is 

We use these expressions to calculate the running coupling constant, i.e. we make no further 
expansion. 

The leading order solution for the valence quark momentum sum rule V2{Q^) using only 
the leading term in the 7-function reads 

W^) = f^V'°vi(Q2). (8) 



«.(«), 



The next to leading order result reads 



^NS f, NS a „]VS 



The NLO solution can be rewritten |jT2|] as the LO solution times a polynomial in as- This 
is the result which we will refer to as truncated NLO. It reads 



^2(Q') = V^iQof as{Q^)/as{Ql 



,(2) 



NS 



(10) 



where the values for d^^ and Zns can be found in Table |. 



For the singlet part the second moment of the distribution function for quarks, S2((5^), 
mixes with the second moment of the gluon distribution, 6*2 (Q^)- One has, however, S2((5^) 
+ G2{Q^) = 1, which makes it possible to write the evolution as 



G2(g')=M^(a,(g^)) 



G2{Ql) 



with 



M-^faJ = exp - / ' da 



, (7gg(«') +lGG{a'y 



(11) 



2[3{a') 



(12) 



Here 7gg(as) and 7gg(q^s) are elements of the singlet anomalous dimension matrix, expanded 
in as in the same way as the nonsinglet anomalous dimension function. For the second 
moment the anomalous dimensions obey = —70^ and 7gg = —IqG- The expansion 
coefficients are 7qg^o = 64/9, 7^^^! = 96.6584 — 10.2716 nj, 7gg,o = 4nj/3 and 7gg,i = 
15.0864 n/. 

The solution for the function M'^(a<j), appearing in the evolution of the singlet quark 
momentum sum rule is in leading order given by 



M^ias 



a. 



2/30 



The NLO solution for M^[as) is given by 



(13) 



M^{as) 



(7qq,0+7GG,0) M-iqq,\+-iGG,l)-fh(-iqq,0 + ^GG,0'> 

as \ "''0 /47r/3o + /5ias(Q^)\ '"""i 



MQD) V4vr/3o + /3i«s(gg) 

This leads to the following LO solution for 6*2, 

(Tg9 ,0+TGG,0) 

G^m = — ^ — + ' 



(14) 



as 



2/3o 



799,0 + 7gg,o \<^s{Qo) 



G2{Ql) 



Iqqfl + iGGfl 



(15) 



For the NLO solution we do not have an analytic expression, but using the result for 

a numerical solution is easily obtained. The truncated NLO result is given by the coupled 

equations 



'i-&)MQl)-&G2{Ql)] 





^(2) 







1 + 



An 



,(2) 



MQl). 

G^iQ') = h(l - a) S2(Qg) + « G^iQl)] 



An 



1 + 



asm - c^siQl 
An 



(16) 





J (2) 







An 



+a-a) 



1 + 



An 



asiQl) 



,(2) 



An 



K 



G 



(17) 



The (i's are the relevant anomalous dimensions for this moment, here evaluated to first order. 
Higher order corrections come from the Z's and i^'s, which are (Q^ independent) coefficients 
tabulated in table |. Note also that conservation of momentum requires S2(Q^)+G'2((5^) = 1, 
which the above moments satisfy due to the relation between a, , and K"^ . 



B. The singlet spin sum rule 

The singlet piece of the first moment of gi has recently been computed to next to lead- 
ing order in [§. This includes both the singlet coefficient function C'^, as well as the 
anomalous dimension, of the singlet axial current in the MS scheme, using dimensional 
regularization. In the adopted normalization, this yields for the sum rule expressed in terms 
of the singlet axial matrix element 

rf(Q') = r gt{x,Q')dx = C'{a.{Q')) AE{Q'), (18) 

J 

with 

"/ 

AS(/i^)s<^ = {p,s\J^\p,s) = {p,s\^qaaj5qi\p,s) = (An + Ad + As H )s„, (19) 

i=l 

the quantity sometimes interpreted as the spin carried on the quarks. The coefficient function 
is given by 
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C^as) = 1 -a,/7r + </7r^ (-4.5833 + 1.16248 n/). (20) 

Not included in the sum rule for gi in Eq. |l^ are higher twist contributions. The matrix 
element in Eq. |18| is scale dependent, 

governed by the anomalous dimension, which with our normalization is 

with 7f = 16 Uf and 7f = 314.67 — 3.556 ri^. 

The LO solution for the singlet axial charge AE reads 



AS(Q2) = exp ( ^(a,(Q2) - aM))) mQl). (23) 



while the NLO solution reads 



fa7f-/307| 



Finally, the truncated NLO solution for AS is 



+Tk^2MQ')-o^sm?]mQl)- (25) 



1287r2/?2 

Experimental results are mostly given for AS(Q^), which is obtained from the exper- 
imental sum rule by explicitly factoring out the coefficient function C'^(as(Q^), but not 
factoring out the dependence in the matrix element, given in Eq. Note that C^{Q'^) 
differs at second order and beyond from the analogous function C^^{Q'^), 

C^^as) = l-as/7r + al/7i^i-4:.5833 + nf/3). (26) 

which appears in the Bjorken sum rule. Being nonsinglet, the Bjorken sum rule has no 
analogous anomalous dimension correction. 



III. RESULTS 



In our previous work |TT|, we proposed fixing a quark model scale, Qq, where e.g. 
G^2(<5o) — 0; and/or V2{Ql) = 1. This can be obtained by evolving from experimental 
values at high |15|. The spin sum rule was considered in the same way, with a starting 
point AE((5o) taken from a quark model value, and then evolved up to relevant to DIS 
experiments. In the bag model the estimate for AE(Qq) ^ 0.65 ||TT[|, the reduction from 
unity coming from the lower components of the (relativistic) quark spinors, the same source 
which reduces the axial charge in the bag model from its nominal value of 5/3. 



The most important improvement of the results of [|lT| is the inclusion of the effects of 
corrections in the next order in a^- This of course does not justify the use of perturbation 
theory in the domain where we are using it, going to rather large values of the strong 
coupling constant, ~ 2. On the other hand, we can get some feeling for the convergence 
or nonconvergence of our approach by comparing first and second order evolution for the 
various moments. 

It is well-known that for the running coupling constant the difference between the first 
and second order results is large when one looks at the functional dependence of on Q^. 
Similarly the evolution of the moments as a function of as can also be strongly dependent 
on the order. Evolving down from as{M'^) = 0.117 P,p!^, and the starting value ^2(4 



GeV^) = 0.44 [T^, using LO order equations, gives G2 = when as = 1.80. Using NLO 
order equations gives G2 = when as = 1.79. Much larger is the difference for the valence 
momentum sum rule. Here one finds that starting from V2(4 GeV^) = 0.40, using LO order 
equations, gives V2 = 1 when as = 2.77. Using NLO order equations gives V2 = 1 when 
as = 2.21. 

When we plot moments against each other as done in Fig |l|, we notice that the NLO 
calculations (dashed line) do not exhibit a drastically different behavior as compared to the 
LO calculations (solid). When we show NLO calculations, the dashed line shows the exact 
solution to the evolution equations. The dotted line shows the NLO truncated expansions 
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in as given in Eqs |T0|, [T^ and |T7[ The dot-dashed hne is the same truncated expansion, 
but any terms involving higher order corrections in are evaluated using the leading order 
expression for as, as suggested in ref [jl2|. Comparison of the dotted, dot-dashed and solid 
line indicate typical uncertainties in the NLO result. The differences between them is one 
order in ag higher. 

The same comparison of LO (solid) and three approximations for the NLO results can 
be seen in Figs ^ and |^, which show AS plotted against V2 and G2 respectively. For the spin 
sum rule, we have compared our results with AS((5cxp = 10 GeV^) because we note from 
Eq. (pID that the results remain proportional to the starting value at all Q^. This makes it 



easy to consider any scenario, e.g. starting from a world average such as AE(4 GeV^) = 0.31 



1^ or starting from a low-energy value [|rT 



In Fig. 0, as is increasing down and to the right. There is no single value of where 
both V2 = 1 and G2 = 0, in either LO or NLO perturbation theory. In both cases G2 vanishes 
earlier (at higher Q^), which is consistent with an intuitive picture that at the quark model 
scale, meson-cloud effects result in some residual qq sea. Note that if one plots e.g. V2 vs as, 
there is a stronger dependence on the order of perturbation theory used, as the evolution of 
as is itself highly modified at these low Q^. It is encouraging that when these observables 
are plotted against one another, the trends are similar. Furthermore, varying the value of 
a5(M|) within current experimental limits (e.g. using values ranging from 0.11 to 0.12) has 
negligible effect on these curves. 

In Fig. 0, we show the relative value of the singlet axial matrix element AS versus V2, 
normalized to the value at Qexp ~ GeV^, the characteristic scale of the SMC experiment. 
The right edge corresponds to a value of Qg where V2{Qq) = 1. We see that AE((5o)/AE(10) 
increases from about 1.71 (LO) to 1.98 (NLO). This increase of the NLO evolution is, unlike 
that in Fig. |1|, sensitive to the value of q;s(M|). Increasing the value of q;s(M|) by 5%, gives 
ratios of 1.83 (LO) and 2.21 (NLO). 

In Fig. 0, AS is plotted against G2, again normalized to its value at Qexp = 10 
GeV^. Evolving from 10 GeV^ to Q^, this time fixed from G2{Qq) = 0, gives for the 
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ratio AS(Q^)/AE(g2^p) values of 1.39 (LO) to 1.68 (NLO). The ratios are smaller than in 
Fig. 0, because G2 = corresponds to a value of V2 < 1. Again a 5% larger value of as(M|) 
results in somewhat larger ratios, 1.45 (LO) and 1.85 (NLO). The persistent enhancement 
of AS at the low-energy scale, also in NLO, suggest that a valence picture with AS of the 
order of 0.5 - 0.8 is consistent with the experimental result in DIS of the order of 0.3 - 0.5. 
Clearly, the scale dependence cannot be neglected in interpreting the results in deep inelastic 
experiments. 

The evolution from = 1 GeV^ to 10 GeV^ is presumably more reliably in the pertur- 
bative regime, and in this case we find a ratio AS(1)/AS(10) of 1.031 (LO), or 1.068 (NLO) 
when as{M^) = 0.117. Especially here, we note a sensitivity to as(M|). Taking its value to 
be 5% higher, the ratios becomes 1.040 (LO) and 1.110 (NLO). The modification of the ac- 
tual singlet moment Ff , which includes C^{Q^) as well, shows a decrease for Ff when going 
to lower momenta, the ratios being Ff(l)/Ff(10) = 0.980 (LO) and 0.956 (NLO). (Taking 
as(M|) 5% higher gives 0.973 (LO) and 0.921 (NLO).) Although fairly small in this region, 
the contribution to the evolution from the anomalous dimension clearly can and should be 
taken into account, and goes beyond the standard QCD effect arising purely from C^iQ"^), 
which is sometimes all that is taken into account. Note furthermore that in Ff((5^), higher 
twist contributions proportional to l/Q"^ could contribute. These have not been included in 
the above estimate for Ff , which refers purely to the twist two part. 

Another point that deserves discussion is the inclusion of the effects of the strangeness 
threshold. If one considers the KK threshold, i.e. ~ 1 GeV^, as an appropriate value, 
a large part of the evolution down to involves Uf = 2. If the strangeness content of 
the nucleon has not become zero, the decoupling of strangeness from the evolution leads to 
ambiguities in the treatment. It would require a global analysis, which takes carefully into 
account existing inequalities such as As{x) < s{x) and the consequences for the moments. 
The evolution equations with nj = 2 instead of = 3 in general tend to somewhat slow 
down the increase of AS at lower momentum scales. 

Finally, if we assume that (i) AS((5o) = 0.65 as an appropriate value at the low mo- 
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mentum scale, e.g. from an effective low energy model |Tl|], and (ii) pQCD continues to 
work down to low Q^, and (iii) the asymptotic values of the nonsinglet combinations of the 
axial matrix elements are known from weak decays, An — Arf = 1.257 and (from low energy 
hyperon decays) Am + Ad — 2As = 0.58 ± 0.05, then we are able to calculate at any scale 
the axial matrix elements for each of the quark flavors. The values at Qq (corresponding to 
6*2 = 0), 1 and 10 GeV^ are given in Table ||. Assuming three active flavors at Qq gives 
a positive value for A^. In this case we have the possibility to incorporate the strangeness 
threshold in a natural way, using only two active flavors to evolve AS down to 0.58, then 
continuing with three active flavors. The numbers in this scenario for NLO are given in Ta- 
ble 0. The strangeness threshold required is Ql = 0.286 GeV^. We note that decoupling of 
the strange quarks in the momentum sumrule at this same threshold implies that at 4 GeV^ 
the momentum carried by the strange antiquarks as compared to nonstrange antiquarks is 
2s2/{u2 + d2) = 0.57. 

In conclusion, using NLO equations, a very satisfactory picture is obtained running all the 
way from a low-energy-scale proton without gluons but with some nonstrange sea, acquiring 
nonzero values for strangeness matrix elements only starting at the strangeness threshold 
which is slightly above the scale where G2 = 0- We have analyzed the errors arising from an 
uncertainty in the strong coupling constant as{Mz) = 0.117 ± 0.005 and those coming from 
an uncertainty in the octet axial charge, 0.58 ± 0.05. These are indicated in Table ||. Note 
that the results for Au and Ad are not sensitive to the octet axial charge if this is taken 
to coincide with the strangeness threshold. Using these results and including second order 
QCD corrections everywhere, we find that 

r?(10 GeV^)= (0.109 ± 0.002) + (0.062 ± 0.007 ± 0.009)AE(g2) 

= 0.149 ±0.006 ±0.006, (27) 

r^(10 GeV^)= (-0.080 ± 0.001) ± (0.062 ± 0.007 ± 0.009)AE(g2) 

= -0.040 ± 0.003 ± 0.006, (28) 

rf(io Gev^)^ o.5(r? ± r^)(i - 1.5 c^^) 
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= (0.013 ± 0.001) + (0.056 ± 0.006 ± 0.008) AE{Ql) 

= 0.049 ± 0.004 ± 0.005 (29) 

(using the usual D-state admixture of 6% in the latter). The first error in each term arises 
from our assumed uncertainties in as{Mz) and from the octet part of the sum rule, here 
added in quadrature. The second error bar (if shown) comes from our estimate of the 
prescription dependence associated with evolving AS from 10 GeV^ to Q^. This includes 
e.g. differences in truncation schemes (see Figs. 2 and 3), choice of s-quark threshold 
mechanism, and determination of Qq from V2 = 1 rather than G2 = 0. These have all been 
discussed above. We conservatively estimate AS((5q)/AS(10) = 1.65 ± 0.25, and the final 
numbers above correspond to this choice, with AS((5o) = 0.65. The relations between the 
experimental sum rules and the "spin carried by the quarks", i.e. AS(Qq) are illustrated in 
Fig. |. 

IV. CONCLUSIONS 

In this paper we have investigated the extent to which measurements of the spin sum 
rule at high energies should be interpreted, in view of the role of their scale dependence. We 
have investigated the spin sum rule together with the momentum sum rules in a systematic 
way, extending our earlier results that used purely leading order evolution to results that 
use next to leading order evolution. This has become possible in part due to the recent 
work of Larin We have estimated uncertainties arising from scheme dependence and 
higher order QCD effects in the evolution, by using several truncation prescriptions. Our 
results indicate that many qualitative features present in the leading order remain the same. 
Quantitative differences show up, but do not upset the picture. In particular, we have 
considered momentum sum rules for valence quarks and gluons compared with one another 
or spin sum rules compared with the momentum sum rules. In general, we have not made 
any interpretation of the relative gluonic contribution to the results for the singlet moment, 
which helps avoid obvious scheme-dependent assumptions, (see, e.g. reference [|l6l, which 
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shows that the gluon contribution is zero in certain renormahzation schemes). We are simply 
using the operator product expansion for the moments of interest, at next to leading order. 

The results for evolution from 1 GeV^ to 10 GeV^ are quite reliable, and the fractional 
change evolving down to Qq is apparently reasonably stable to next to leading order cor- 
rections too. Of course, the moderate sensitivity in going from first to second order is no 
proof of the reliability of the perturbative expansion. We have noted that while the sensi- 
tivity to using first or second order is moderate, the sensitivity to the value of as(M^) is 
quite strong. We find that the inclusion of a strangeness threshold is not very important 
for the rate of evolution. Much more important is the fact that there is a large change in 
AS running from 10 GeV^ to a low-momentum scale Qq, large enough that it can easily 
reach a point where AS = Am + Ad — 2 As, i.e. As = 0, a natural point for the strangeness 
threshold. We have made this more explicit by starting with a value |jll| of AT,{Ql) = 0.65, 
although any assumption that one can estimate AS(Qq) from a low energy nucleon model 
is very clearly subject to debate. This scenario implies that at present there is no reason to 
require an anomalously large strange contribution in the proton at a low-momentum scale. 
Of course, independent measurements of the strangeness content at this scale (e.g. from 
elastic neutrino - nucleon scattering) are important to confirm this. 
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FIGURES 

FIG. 1. Plot of G2 vs V2. Solid (dashed) curve shows first(second) order calculations. We 
started from experimental values ||l5| at roughly 4 GeV^ and evolve downwards. Dotted curve is 
2nd order, but truncated. Dash-dotted is again 2nd order, truncated, but with as in all higher 
order terms replaced with its leading order expression. 

FIG. 2. Same as Fig. 1, but plotting AS/AS(Qg) versus V2. Curves are labeled as before. 

FIG. 3. Same as Fig. 3, but plotting AS/AS(Qg) versus G2. 

FIG. 4. Experimental sum rules as a function of the spin carried by the quarks, A$]((5q), 
including uncertainties from a^, octet axial charge, and scheme dependence are given by the shaded 
areas. The areas enclosed by the dashed lines are the results as a function of AS(10 GeV^). 
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TABLES 

TABLE L Numerical values [p!^Jl^ used for the various parameters appearing in the truncated 
next to leading order solutions for the momentum sum rules, Eqs 10 through 



Uf 


a 


"NS 


4'^ 


z+ 








2 


0.2727 


0.3678 


0.5058 


1.486 


1.428 


0.4544 


-0.1704 


3 


0.36 


0.3951 


0.6173 


1.783 


1.507 


2.121 


-1.193 


4 


0.4286 


0.4267 


0.7467 


2.355 


1.654 


5.895 


-4.421 


5 


0.4839 


0.4638 


0.8986 


3.341 


1.904 


22.604 


-21.191 



TABLE IL Values for the axial matrix elements of the quarks using a starting value of AS((5o) 
= 0.65 at the scale where G2{Qq) = using NLO results and a strangeness threshold at the point 
where As = 0. The errors refer to uncertainties in as{Mz) = 0.117 it 0.005, and in the octet axial 
charge, 0.58 it 0.05 (underlined errors). 



Q2 [GeV^] Au Ad As AS 

Q§ 0.954 -0.304 - 0.650 

q2 ^ 0.3 =F 0.1 0.919 ± 0.025 -0.339 ± 0.025 0.00 0.580 ± (LOS 

1 0.873 T 0.01 -0.384 ^ 0.01 -0.045 T 0.01 =F 0.025 0.445 ^ 0.03 T 0.025 

4 0.866 T 0.01 -0.391 ^ 0.01 -0.052 T 0.01 =F 0.025 0.423 ^ 0.03 T 0.025 

10 0.864 T 0.01 -0.393 T 0.01 -0.055 T 0.01 T 0.025 0.417 T 0.03 T 0.025 
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This figure "figl-l.png" is available in "png" format from: 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9412287vl 
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This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9412287vl 
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This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9412287vl 
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